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Abstract 

The perpendicular and parallel multilayer relaxations of fee (210) surfaces are studied using equivalent crystal 
theory (ECT). A comparison with experimental and theoretical results is made for AK210). The effect of 
uncertainties in the input parameters on the magnitudes and ordering of surface relaxations for this semiempirical 
method is estimated. A new measure of surface roughness is proposed. Predictions for the multilayer relaxations and 
surface energies of the (210) face of Cu and Ni are also included. 


1. Introduction 

In the last few years there has been consider- 
able interest in the study of the surface structure 
of high-index faces of metals [1-10]. Experimen- 
tal evidence of both perpendicular (to the surface 
plane) and parallel (without loss of symmetry) 
relaxations in several open metal faces, provided 
the necessary background for theoretical studies 
[ 1 1 — 16]. Following the calculations of Barnett et 
al. [13], which first predicted the occurrence of 
significant interlayer spacing and registry relax- 
ations, the surface structure of six fee Al surfaces, 
including the (210) face, was studied by Jiang et 
al. [11] using the modified point-ion model. A 
low-energy electron diffraction (LEED) analysis 


* Corresponding author. 


by Adams et al. [3] provided the only experimen- 
tal evidence to date of perpendicular and parallel 
relaxation of Al(210) surfaces. More recently, 
Sinnott et al. [10] reported corrected effective 
medium (CEM) results for the surface energies of 
several fee metals and provided theoretical esti- 
mates for the perpendicular relaxation of the first 
two interplanar spacings. Finally, Chen and Voter 
[12] performed a calculation using embedded- 
atom potentials raising the issue of the possibility 
of reconstruction in fee (nlO) (n = 1, 2, 3) sur- 
faces although none was found experimentally in 
such cases [3,9]. These theoretical studies, to- 
gether with several other first-principles or 
semiempirical calculations for the determination 
of metallic surface structure, have provided not 
only large amounts of data but also some insight 
on the different mechanisms involved. However, 
some issues which could be considered somewhat 
minor given the success of most of these studies, 
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have not enjoyed the same level of attention, the 
focus being in all cases up to what extent these 
techniques were able to reproduce the experi- 
mental results. 

In this paper, we intend to open the discussion 
on some of these issues, by focusing on a multi- 
layer relaxation calculation of the Al(210) surface 
by means of a semiempirical method which, in 
general, has been very successful in previous ap- 
plications to surface phenomena, namely, equiva- 
lent crystal theory [14-16]. As opposed to the 
case of low-index surfaces, where ECT provides, 
in general, good results, the predictions for the 
AK210) surface are not as good. Therefore, it 
serves as a useful example for answering the 
question of what ingredients of the theory are 
relevant, and which ones are not, in order to 
provide a physically accurate description of the 
problem at hand. 

Like any of the methods available for the 
theoretical study of surface phenomena, whether 
it is a first-principles calculation or a semiempiri- 
cal one, approximations are made and some ex- 
ternal input is used, without a clear understand- 
ing of how these assumptions translate into the 
final results. Moreover, semiempirical techniques, 
which became the standard for simulations, over- 
simplify the problem in the attempt of making the 
computational aspects simple enough to allow for 
lengthy and complex calculations. Often, it is at 
the expense of introducing mechanisms whose 
roles, partly because of their generality, become 
obscure in that their influence on the final results 
is not well understood. The criteria used to inter- 
pret the results is also important: while the out- 
put may consist of a single quantity, easily identi- 
fiable (in the case of surface relaxation, the per- 
centage changes in interplanar spacing), it is of- 
ten the case that alternative solutions can still 
shed some light on the behavior of the system, 
even if they do not correspond to observed exper- 
imental results (i.e., alternative relaxation pat- 
terns which could correspond to local minima in 
the multidimensional energy surface). We there- 
fore focus our attention on three aspects of the 
surface structure calculation: the input data used 
and its influence on the results, the energy algo- 
rithm and its ability to deal with the main physi- 


cal effects that take place in the system, and the 
context in which the output data is analyzed. 

In the particular case of ECT, although some 
of these features are common to other semiem- 
pirical techniques, it is important to analyze how 
uncertainties in the experimental data used as 
input affect the outcome of the calculation and, 
in the problem of multilayer relaxation, how that 
influences the ensuing relaxation pattern. We will 
find that very small uncertainties in the input 
data used in ECT (which is also the input of 
other semiempirical techniques) are amplified 
into large uncertainties in the final results. Sec- 
ondly, we concentrate on the role of each of the 
terms that enter in the calculation of the energy 
by studying the results obtained under different 
parameterizations and therefore gaining some in- 
sight into their physical interpretation. In the 
study of multilayer relaxation, this relates to the 
mechanisms provided by the energy algorithm to 
describe the bond length anisotropies which ulti- 
mately dictate the structure of the surface. Fi- 
nally, we generalize the concept of surface rough- 
ness thus establishing a more appropriate frame- 
work for the analysis of surface relaxation pat- 
terns theoretically obtained. 

We organize the paper as follows: in Section 2 
we briefly discuss the concepts of equivalent crys- 
tal theory. In Section 3 we concentrate on the 
first of the three main issues of this paper by 
introducing theoretical ‘error bars’ to the predic- 
tions of low-index surface relaxation, which we 
later generalize in Section 4 to the case of high- 
index surfaces. We discuss the dependence of the 
ECT results on the parameterization chosen for 
the bond-compression term in the ECT energy 
expansion. We also introduce a new concept to 
replace the definition of surface roughness and 
analyze the ECT results in this new framework. 


2. Equivalent crystal theory 

Equivalent crystal theory [15,16] is based on an 
exact relationship between the total energy and 
atomic locations and applies to surfaces and de- 
fects in both simple and transition metals as well 
as in covalent solids. Lattice defects and surface 
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energies are determined via perturbation theory 
on a fictitious, equivalent single crystal whose 
lattice constant is chosen to minimize the pertur- 
bation. The energy of the equivalent crystal, as a 
function of its lattice constant is given by a uni- 
versal binding energy relation [17]. 

Let € be the total energy to form the defect or 
surface, then 

e=I>,, (1) 

/ 

where e, is the contribution from an atom i close 
to the defect or surface. ECT is based on the 
concept that there exists, for each atom /, a 
certain perfect, equivalent crystal with its lattice 
parameter fixed at a value so that the energy of 
atom / in the equivalent crystal is also c f -. This 
equivalent crystal differs from the actual ground- 
state crystal only in that its lattice constant may 
be different from the ground-state value. We 
compute € ( via perturbation theory, where the 
perturbation arises from the difference in the ion 
core electronic potentials of the actual defect 
solid and those of the effective bulk single crystal. 

For the sake of simplicity, the formal perturba- 
tion series is approximated by simple, analytic 
forms which contain a few parameters, which can 
be calculated from experimental results or first- 
principles calculations. Our simplified perturba- 
tion series for e, is of the form 

e, = AE{F*[a*(0] + £*"[«!(/,;)] 

' j 

j.k 

+ ( 2 ) 

pa ' 


where 


F*[a*] = 1 -(1 +a*)e" fl * (3) 

and A E is the cohesive energy. In this expres- 
sion, we distinguish four different contributions 
to the energy of atom i and thus, the existence of 
four different equivalent crystals which have to 
be determined for each atom /'. 

The first term, F*[a*(i)], contributes when 
average neighbor distances are altered via defect 
or surface formation. It can be thought of as 
representing local atom density changes. In most 
cases, this ‘volume’ term is the leading contribu- 
tion to c, and in the case of isotropic volume 
deformations, it gives e J to the accuracy of the 
universal energy relation [17]. The value of a*(i\ 
the lattice parameter of the first equivalent crys- 
tal associated with atom /, is chosen so that the 
perturbation (the difference in potentials be- 
tween the solid containing the defect and its bulk, 
ground-state equivalent crystal) vanishes. Within 
the framework of ECT, this requirement trans- 
lates into the following condition from which 
af(0 is determined: 


NR f exp( — a/?, ) + MR% exp 



*2 


- £ rf exp{— [« + S(r ; )]r y } =0, (4) 

defect 


where the sum over the defect crystal or surface 
is over all neighbors within second-neighbor 
(NNN) distance. r J is the actual distance between 
atom i and a neighbor atom /, N and M are the 
number of nearest-neighbors (NN) and next- 
nearest-neighbors, respectively, of the equivalent 
crystal (12 and 6 for fee, 8 and 6 for bcc) and /?, a 
and A are parameters known for each atomic 


Table 1 

Computed constant and experimental input for ECT 


Element 

P 

/ 

a 

A 

10 ~ 2 A 7t /D 

Ur'AJD 

10 4 d 

AE 

fl c 

Al 

4 

0.336 

2.105 

0.944 

7.822 

2.104 

591.4 

3.34 

4.05 

Cu 

6 

0.272 

2.935 

0.765 

5.784 

2.530 

99.74 

3.50 

3.615 

Ni 

6 

0.270 

3.015 

0.759 

7.382 

2.793 

100.1 

4.435 

3.524 

Fe 

6 

0.277 

3.124 

0.770 

9.183 

1.887 

60.62 

4.29 

2.86 


The constant p is 2n - 2, where n is the atomic principal quantum number, / (in A) is a scaling length and A (in A) is a screening 
parameter (see text) The constants A 3 and A 4 are dimensionless. AE (in eV) is the cohesive energy and a c (in A) the equilibrium 
lattice constant. 
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species. Table 1 displays the values of these pa- 
rameters for the fee elements used in this work 
(see Ref. [16] for a complete list). S(rj) is a 
screening function and and R 2 are the NN 
and NNN distances in the equivalent crystal. The 
equivalent lattice parameter, a l9 is thus related to 
the scaled quantity a* via 



where r WSE is the equilibrium Wigner-Seitz ra- 
dius, / is a scaling length and c is the ratio 
between the equilibrium lattice constant and 

r WSE* 

The higher-order terms are relevant for the 
case of anisotropic deformations. The linear inde- 
pendence attributed to these four terms is consis- 
tent with the limit of small perturbations which 
we assume for the formulation of ECT. The sec- 
ond term, F*[a*(i, ;)], is a two-body term which 
accounts for the increase in energy when NN 
bonds are compressed below their equilibrium 
value. This effect is also modeled with an equiva- 
lent crystal, whose lattice parameter is obtained 
by solving a perturbation equation given by 

NR? exp( -aRj) - NRfi exp( -aR 0 ) 

+ A 2 R$ E (Rj - R 0 ) exp [ -P (Rj - /?„)] = 0, 

( 6 ) 

where (3 = 4 a for the metals used in this work, 
and Rj is the NN distance of the equivalent 
crystal associated with the deviation of NN bond 
length Rj from R 0 , and R 0 is the bulk NN 
distance at whatever pressure the solid is main- 
tained (generally, R 0 is the ground-state, zero- 
pressure value). A 2 is a constant determined for 
each metal (see Table 1 for a list of values of A 2 
used in this work). The scaled equivalent lattice 
parameter is then 

fl 2 = |“ ^*WSE | /^* (7) 

The third term, j, k)] accounts for 

the increase in energy that arises when bond 
angles deviate from their equilibrium values of 
the undistorted single crystal. This is a three-body 


term and the equivalent lattice parameter associ- 
ated with this effect is obtained from the pertur- 
bation equation 

NR f exp( -a/?,) - NRfi exp (~aR 0 ) 

+ A 3 Rg exp[ -a(R ] +R k ~ 2/?„)] 

sin(0,* - 0) = 0, (8) 

where A 3 is a constant listed in Table 1 and 0 jk is 
the angle between the NN distances Rj and R k 
with the atom i at the center. 0 is the equilibrium 
angle, 70.5° for bcc and 90° for fee. This term 
contributes only when there is a bond-angle 
anisotropy ( 6 jk 0). The scaled lattice parameter 
is then 

*!-(" ^WSE | /l ’ (9) 

The fourth term, F*[a*0‘, p, q)\ describes 
face diagonal anisotropies (see Ref. [16] for a 
detailed description, for each lattice type, of the 
structural effect associated with this term). The 
perturbation equation reads 

NR? exp( —aR x ) - NR$ exp (~aR 0 ) 


exp[-a(Rj + R k + R, + R m - 4R 0 )\ =0, 

( 10 ) 

where d is the face diagonal of the undistorted 
cube and A 4 is a constant adjusted to reproduce 
the experimental shear elastic constants (Table 
1). Finally, 

a * = ( ” r WSE j/^- 

Consider a rigid surface (i.e., no interlayer 
relaxation): all bond lengths and angles retain 
their bulk equilibrium values, thus F*(a*) = 
F*(a*) = F*(a*) = 0. The surface energy is 
therefore obtained by solving for the ‘volume’ 
term represented by F*(af ) only. If we consider 
a rigid displacement of the surface layer towards 
the bulk, as is the case in most metallic surfaces, 
the higher-order terms become finite: some bonds 
are compressed, contributing to F*(^*), the bond 
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angles near the surface are distorted as well as 
the difference between face diagonals in some 
cases, generating an increase of energy via F*(a| ) 
and F*(a *). For the case studied in this work, 
these additional contributions to e t are generally 
small, representing only 1% to 2% of the total 
energy. However, while these anisotropy terms 
are small for metals when there is no reconstruc- 
tion, they play an important role in the energetics 
of these defects where the differences in energy 
between the rigid and relaxed configurations are 
also small. In what follows, we will refer to this 
ECT formalism as ECT II. 

An earlier version of ECT [15], which we will 
refer to as ECT I, provides a simpler, although 
less accurate framework for a defect calculation. 
The second term in Eq. (2) is replaced by a 
simple expression, which allows for the direct 
calculation of the energy associated with bond- 
compression effects, 

^ M n 6 

e 2 = A£ E E j^F*(a* mn ), (11) 

where N s is the number of atoms in the solid, 
0 mn = 1 if n < 0 and 8 mn = 0 otherwise, M n is 
the number of nearest neighbors of atom n , L mn 
is the number of nearest neighbors of atom m or 
n, whichever number is smaller, and a* m n is given 
by 


a 


* 

mn 


R m n/ C l r WSE 

/ 


( 12 ) 


with 


/ = 


A E 

\2tt Bryy 


(13) 


B is the bulk modulus of the crystal, R mn is the 
distance between atoms m and n , c, is the ratio 
of the equilibrium nearest-neighbor distance in 
the crystal to r WSE , and r WSE is the equilibrium 
Wigner-Seitz radius. In ECT I [15], the third and 
fourth term of the energy expansion (Eq. (1)) are 
ignored. In what follows, we will list results as 
obtained with either one of the two versions of 
ECT. Those results obtained with the full energy 
expansion (ECT II) [16], will be analyzed in terms 
of the value of the parameter (3 which dictates 


the ‘strength’ of the bond-compression term 
therefore playing an important role in the ener- 
getics of surface relaxation as it will be seen that 
this term is mainly responsible for avoiding the 
collapse of the top layers onto each other. 


3. Uncertainties in the prediction of multilayer 
relaxation 

Before proceeding to the calculation of multi- 
layer relaxation in high-index faces, we will dis- 
cuss some features of theoretical calculations of 
these quantities. Ref. [14] provides a reasonably 
large sample of both experimental and theoretical 
results for changes in interlayer spacing in pure 
fee and bcc crystals. In all cases, the semiempiri- 
cal, theoretical techniques used, rely either on 
input data (generally experimentally determined) 
or on certain approximations for some of the 
variables of relevance. Necessarily, results will 
depend on such choices. Multilayer relaxations 
involve at best very small changes in position, 
with correspondingly, small changes in surface 
energy, whose minimization is the criterion used 
to determine the final interlayer spacings. Thus, 
the search for a minimum of the surface energy, 
as accurate as the minimization technique might 
be, will be strongly influenced by the approxima- 
tions made, the error in input parameters and the 
shallowness of the minimum in the surface energy 
surface. As a consequence, to quote just one 
value for each of the changes in interlayer spac- 
ings as is ordinarily done, might not reflect the 
ambiguities in these calculations. In this paper we 
adopt a different path: to each theoretical predic- 
tion, we will attach an estimate of the possible 
errors due to any of the reasons mentioned above. 
Although there is no certain way to determine 
such errors (after all, the predictions are, within 
their own framework, exact), we will see that 
changes on the order of 1% in the surface energy 
can generate quite interesting variations in the 
relaxation schemes predicted. In particular, within 
the framework of ECT (I [15] or II [16]), such 
small changes in the surface energy can be easily 
obtained by changing any of the input parameters 
(lattice constant, cohesive energy, bulk modulus) 
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by a similar amount, well below the usual experi- 
mental errors in the determination of such quan- 
tities. 

To illustrate this issue, we will focus our atten- 
tion on the surface structure of some fee pure 
metals (Al, Cu and Ni). As can be seen in Tables 
2-11 of Ref. [14], previous theoretical and experi- 
mental studies show a wide spread in the predic- 
tions of the changes in interlayer spacings for the 
(100) and (110) surfaces. Even results obtained 
within the same theoretical technique (em- 
bedded-atom method (EAM) [22], ECT) do not 
agree with each other (due to different fitting 
procedures of the embedding function in the case 
of EAM and different input data in both cases). 
Although there is general qualitative agreement, 
regarding the contraction or expansion pattern 
found for successive layers, in some cases the 
absolute theoretical values predict the wrong 
trend with respect to experimental results (see, 
for example, A1(1()0)). The ECT II results (from 
Refs. [14] and [16]) also highlight this inconsis- 
tency. The difference between the values ob- 
tained in this work and those from previous ap- 
plications of ECT is easily traceable to slightly 
different values of some of the input parameters. 

As mentioned above, in order to account for 
these and other ambiguities in the calculation, we 
investigated the change in predicted relaxations 
due to small changes in the rigid surface energy. 
We thus defined ‘error bars’ in such way that all 
the intermediate values so obtained predict varia- 
tions in surface energies within a certain toler- 
ance. In this work, we set the tolerance at 1% of 
the equilibrium surface energy cr e . This defines a 
surface cr(A d 12 , Ad 23 ) and the allowed values for 
these parameters are such that <r e < cr(A<i 12 , 
Ad 23 ) < 1.01cr c . Needless to say, this range of 
values does not include all the possible sets (A */ 12 , 
A d 23 ) that correspond to surface energies within 
the allowed values. It is interesting to note, how- 
ever, that in most cases, all the experimental as 
well as theoretical predictions fall within the range 
of uncertainties in such procedure. This defini- 
tion of the error bars is, of course, arbitrary and 
it was chosen as a means to simply illustrate the 
influence of uncertainties in the final results. 

It should be noted that when comparing our 


theoretical predictions with available experimen- 
tal results, the error bars quoted from experiment 
or theory are similar in that the optimum relax- 
ations are determined by minimization of some 
property by varying the input parameters. To 
illustrate this point, we first discuss the surface 
energies and multilayer relaxations of the unre- 
constructed low-index surfaces of pure Al, Ni, 
and Cu crystals. In Table 2 we display the ECT II 
predictions for the surface energies and compare 
the results with typical experimental values for 
polycrystalline samples [18,19]. We note that ex- 
perimental values for the surface energies are for 
polycrystalline surfaces, thus could be strongly 
dominated by the predominant surface plane. The 
experimental values from Ref. [21] have the ad- 
vantages of the data being taken on solids (in- 
cluding low-temperature values), and the data 
being in much better agreement with modern, 
first-principles calculational results. 

In Table 3 we compare results for the multi- 
layer relaxations of the first two interlayer spac- 
ings for those cases for which recent experimental 
data are available [23,31]. The inclusion of the 
theoretical ‘error bar’, as mentioned above, al- 
lows for a better comparison with experiment as 
it shows that for most cases, small changes in the 
input parameters of the method may account for 
the whole range of possible experimental results. 
The exceptions are Al(lOO) and Al(lll), where 
the outward relaxation of the surface layer has 
been attributed to an electron promotion effect 
[32], Semiempirical methods (ECT, EAM, etc.), 
unless specifically designed to do so, do not gen- 
erally allow for such fine electronic structure 


Table 2 

Experimental (Exp.) and relaxed ECT II surface energies of 
Al, Cu and Ni (in erg/cm 2 ) 


Technique 

Al 

Cu 

Ni 

Exp. [18] 

1200 

1790 

2270 

Exp. [19] 

1140 

1780 

2380 

Exp. [19] 

1180 

1770 

2240 

Exp. [21] 

1169 

2016 

2664 

ECTU00) 

1203 

2309 

2982 

ECTU10) 

1284 

2373 

3073 

ECT(lll) 

856 

1767 

2274 



210 
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Table 3 


Surface relaxations of Al, Cu and Ni as percentages of the bulk interplanar spacings 


Element 

Face 

Experiment 



ECT II 

ECT II (two-layers) 


\d n 

A d n 

Ref. 

&d l2 

Arf, 2 

a</ 23 

Al 

(100) 

+ 1.8 


[22] 

-4.68 ± 1.62 

-5.05 ± 1.58 

+ 3.35 ±0.80 


(110) 

-8.5 ± 1.0 

+ 5.5 ± 1.1 

[23] 

-8.29 + 2.35 

-9.53 ±3.58 

+ 1.90 ±2.24 


(111) 

+ 1.7 ±0.3 

+ 0.5±0.7 

[24] 

-3.67 ± 1.21 

- 3.94 ± 1.19 

+ 2.75 ±0.61 

Ni 

(100) 

-3.2 + 0.5 


[25] 

-3.53 + 1.68 

-3.82 ± 1.68 

+ 2.48 ± 0.85 


(110) 

— 9.0 ± 1.0 

+ 3.5 ± 1.5 

[26] 

-6.32 ±2.44 

-6.55 ±3.63 

+ 0.34±2.24 


(111) 

-1.2+ 1.2 


[27] 

-2.89 ± 1.29 

-3.10 ± 1.25 

+ 2.12 ±0.63 

Cu 

(100) 

-2.1 

+ 0.45 

[28] 

-3.52+ 1.74 

-3.81 ± 1.70 

+ 2.47 ± 0.86 


(110) 

-7.5 + 1.5 

+ 2.5 ± 1.5 

[29] 

-6.31 ±2.46 

-6.51 ±3.83 

+ 0.29 ±2.44 


(111) 

-0.7 + 0.5 


[30] 

-2.88 ± 1.30 

-3.10 ± 1.25 

+ 2.12±0.63 


The ECT II &d ]2 column displays results for relaxations of the top layer only while the ECT II (two layers) columns display results 
for the case when the top two layers are allowed to relax. 


effects, thus it is not surprising that our results 
for A d x2 in these cases predict surface layer con- 
tractions, even when the 'error bar is taken into 
account. For completeness we also include results 
for the surface relaxation when only the top plane 
is allowed to relax, in order to single out correla- 
tions with subsequent interlayer spacing changes 
on the surface plane. Again, the agreement with 
available experimental data is good in all cases. 


4. Results and discussion 

As mentioned above, the bond length 
anisotropy term is of special importance as it 
models, in a straightforward fashion, the stiffness 
of the 'springs’ between atoms, thus dictating the 
final configuration of the surface. The bond-com- 
pression term, as defined in Eq. (29) of Ref. [16], 
has a material-dependent parameter ((3) which, 
although fixed in the original formulation, can be 
freely varied. In this work, in addition to using 
the original value of (3, we studied the depen- 
dence of the interlayer relaxations and registry 
shifts on this parameter (3. We also quote results 
obtained with ECT I, for a comparison of its 
effect on relaxation. Finally, being that this work 
is a continuation of our previous work on multi- 
layer relaxation of high-index fee and bee surfaces 
[14], we will also compare our results with those 
quoted in Ref. [14], in order to highlight the 


effect of the inclusion of parallel relaxations in 
the energetics of fee (210) surfaces. 

We will now focus our attention on the Al(210) 
surface. Table 4 shows the results for the surface 
energy of Al(210) obtained with different ap- 
proaches: ECT I [15], ECT II [16], embedded- 
atom method [22] and corrected effective-medium 
theory [10]. 

The different entries for ECT correspond to 
the unrelaxed case (“Rigid”), the perpendicularly 
relaxed case as obtained with the earlier version 
of the bond-compression term [15] (“ECT I ± ”) 
and the current (ECT II) version [16] (“ECT II-4 
1 ”) for the stiffness parameter (3 = 4 a, with a 
defined in Ref. [20], and the case with perpendic- 


Table 4 

Surface energy of the AJ(210) face (in erg/cm 2 ) (see text) 


Method 

_ AJ 
a 210 

Rigid 

1493 

ECT I X 

1426 

ECTII-4 X 

1405 

ECT I x + II 

1424 

ECT 11-4 X + || 

1404 

ECT II-3 X + || 

1390 

ECT II-2 X + I! 

1369 

CEM-LMTO 

1440 

CEM-EMP 

1330 

MD/MC-CEM-LMTO 

1150 

MD/MC-CEM-EMP 

1110 

EAM 

999 
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ular and parallel relaxation using the ECT I 
bond-compression term (“ECT I J_ + ||”) and 
the current version (“ECT II-/? _L + II”) for (3 = 
na ( n — 2, 3, 4). The CEM entries are labeled 
according to the embedding functions used: 
semiempirical (“CEM-EMP”) and those obtained 
from linear-muffin-tin-orbital calculations 
(“CEM-LMTO”), which is their best estimate. 
We quote results from Ref. [10] where the au- 
thors use approximate versions of CEM, used in 
molecular dynamics and Monte Carlo simulations 
(“ MD /MC-CEM-LMTO” and “MD/MC-CEM- 
EMP”) again with two different choices for the 
embedding functions. Substantial differences ex- 
ist between the various results for the surface 
energies. We also include embedded-atom results 
as obtained by Chen and Voter [12]. These results 
are compared to experimental values of polycrys- 
talline A1 samples [18-21], which should corre- 
spond to an average of its highest density planes. 

Table 5 reproduces the results obtained with 
the different variations of ECT described above, 
indicating the contributions from the different 
many-body terms included in surface energy a. ^ 
(i = 1, 2, 3, 4) indicate the /-body contribution to 
the energy. As is the case for perpendicular relax- 
ations [14,16], the absence of reconstruction ef- 
fects translates into <r 3 and a A being much smaller 
than the first two contributions to the surface 
energy a. This is an important issue, as a compar- 
ison between ECT I and II is justified only when 
cr 3 and a 4 are very small. 

When comparing the results for relaxed sur- 
faces, we distinguish between those calculations 
that include parallel relaxations and those that do 
not. Table 6 displays results for the unrelaxed 


Table 5 

Surface energy of AK210) (<r) (in erg/cm 2 ) and the different 
ECT contributions (see text) 


ECT 

a 

0"i 

a 2 


<?4 

ECT II-4 

1404.31 

1369.44 

32.57 

1.40 

0.90 

ECT II-3 

1390.46 

1348.24 

39.12 

1.75 

1.35 

ECT II-2 

1369.33 

1315.54 

49.25 

2.19 

2.36 

ECT I _L 

1426.36 

1371.13 

53.26 

0.71 

1.26 

ECT I ± + 1! 

1424.25 

1366.83 

54.90 

1.55 

0.89 


case and when only perpendicular relaxations of 
the interplanar spacings are included. Table 7 
concentrates on the fully-relaxed case, for which 
LEED experimental data are available [3]. We 
also include results of pseudopotential calcula- 
tions by Barnett et al. [13], and the values for the 
forces on the surface layers of the unrelaxed 
structure made using the point-ion model com- 
puted by Adams et al. [3], which can be taken as 
a representation of the trends of the relaxations. 
In this last set of results, it was assumed that the 
relaxations are linearly proportional to the forces 
in the limit of small relaxations and the actual 
forces are multiplied by an arbitrary factor in 
order to obtain numerical agreement with the 
value of Ad l2 . 

The uncertainties in the experimental values, 
which are the results of a multivariable fit to 
LEED theory, for the parallel relaxations are 
large enough to make it difficult to extract a 
relaxation pattern to which theoretical predic- 
tions can be compared. If we are to take the 
trends of these results seriously, ECT predicts 
different trends from the experimental values for 
Al(210). However, certain features are common 


Table 6 

Percentage change in interlayer spacing perpendicular to the surface of Al(210) 


1 

ECT 


CEM 


MD/MC-CEM 

ECT I _L 

ECT 11 4 _!_ 

LMTO 

EMP 

LMTO 

EMP 

A d l2 

-9.5 

-8.2 

-19.7 

-13.8 

13.8 

-18.9 

Arfy 

-5.6 

-7.5 

-1.0 

+ 3.2 

-1.0 

+ 2.8 

A<A.4 

4-0.8 

+ 2.1 





A<*4.5 

-3.3 

-3.7 





Arf 5 .6 

+ 2.7 

+ 4.0 
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Table 7 

Percentage change in interlayer spacing (A d tJ ) and registry shift (A a (J ) of AK210) 


±11 

Experiment 
Ref. [3] 

Barnett et al. 
Ref, [13] 

ECT 




P1M 

Ref. [3] 

ECT I 

ECT II-4 

ECT II 3 

ECT II-2 

A 

-15.5 ±2.4 

-27.7 

-10.4 

-8.1 ±4.4 

- 10.9 

-13.6 

-15.5 

A ^2,3 

-0.8 ±2.9 

-10.2 

-5.8 

-7.1 ±3.8 

-8.2 

-10.6 

-3.0 

Arf 3,4 

+ 8.9 ±2.6 

+ 25.9 

+ 1.4 

+ 2.9 ±4.2 

+ 2.4 

+ 2.9 

+ 1.5 

^4,5 

-4.4 ±3.6 

-12.8 

-3.9 

-3.4 ±5.8 

-4.1 

-3.5 

-0.3 

^5,6 

-1.2 ±4.6 

-2.4 

+ 3.0 

+ 4.2 ±6.8 

+ 4.5 

+ 4.7 

+ 0.1 

Ao 12 

-0.1 ±3.4 

-2.5 

-0.3 

-0.2 ±2.4 

-0.4 

-0.4 

+ 1.9 

A a 2,3 

-3.2 ±3.1 

-10.0 

+ 0.1 

+ 0.0 ±2.6 

+ 0.1 

-0.1 

-3.2 

A«3.4 

+ 1.7 ± 3.1 

+ 3.8 

+ 1.2 

+ 0.8 ±2.9 

+ 0.9 

+ 0.9 

+ 1.9 

A«4,5 

-2.0 ±4.0 

-1.0 

+ 0.1 

+ 0.0 ±3.5 

0.0 

+ 0.1 

-0.4 

A "5.6 

-0.9 ±4.6 

-2.0 

-0.6 

-0.5 ± 4.4 

-0.7 

-0.6 

0.0 


to all the results: a contraction of d n and d 23 , an 
expansion for d 34 and a contraction of d 45 . The 
theoretical predictions for d 5f) indicate a small 
expansion, which is not in complete disagreement 
with the experimental result, given the large un- 
certainty quoted. Although the trends agree, the 
ECT II-4 results seem to predict the magnitude 
of the contractions incorrectly. We can see (Table 
7) that modifying this term allows the possibility 
of improving quantitative agreement with experi- 
ment. Such modifications await further experi- 
mental results with smaller scatter. It would be 
interesting to test the theoretical results for ECT 
and other theoretical predictions for the corre- 
sponding LEED /^-factors in order to make a 
direct comparison with experimental results. The 
difference in the magnitude of the relaxations 
goes beyond the simple numerical appearance 
as they translate into quite different atomic re- 
arrangements: while the experimental values sug- 
gest a ‘filling’ of the space between surface atoms, 
the ECT solution describes a highly symmetric 
distribution where the same effect is obtained by 


Table 8 

Rorocity of low-index unrelaxed fee faces 

( 100 ) ( 110 ) ( 111 ) ( 210 ) 


B, 

1.2732 

1.8006 

1.1026 

2.8467 

b 2 

1.0922 

1.1691 

1.0426 

1.6116 

b 3 

1.0922 

1.1691 

1.0426 

1.2820 

B, 

1.0922 

1.1691 

1.0426 

1.2360 


a larger net ‘motion’ of the surface atoms toward 
the bulk. 

At this point, we find it convenient to extend 
the concept of roughness of a surface [2] by 
defining the borocity of a surface as 

A t 

B p = ~p , ( 14 ) 

E + Q~ z ' /a 

i- 1 

where A T is the unit surface area and A { is the 
fraction of A r including the projections of the 
hard spheres representing the atoms in layer /, of 
radius half the nearest-neighbor distance, not 
covered by similar ‘disks’ in layers above layer i. 
z i is the location of layer i measured perpendicu- 
larly to the surface plane and a is the equilibrium 
lattice parameter (so defined, B { is the usual 
roughness of the surface [2]). This quantity pro- 
vides a better measurement of the electronic 
smoothness of a surface by attempting to include, 
in a simple fashion, the contributions to the sur- 
face electron density by atoms below the surface 
plane. If the borocity of a surface is to be taken 


Table 9 

Borocity of the relaxed AK210) face 



Rigid 

Experiment 

ECT I 

ECT 11 

B t 

2.8467 

2.8467 

2.8467 

2.8467 

b 2 

1.6116 

1.5873 

1.5984 

1.5957 

B, 

1.2820 

1.2643 

1.2634 

1.2619 

B* 

1.2360 

1.2149 

1.2177 

1.2159 



G. Bozzolo et al. / Surface Science 315 (1994) 204-214 


213 


Table 10 

Surface energy of the Cu(210) and Ni(210) face (in erg/cm 2 ) 
(see text) 


Method 

02io 

(Y Ni 
°210 

Rigid 

2618 

3407 

ECT I X + II 

2559 

3330 

ECT II-4 X + || 

2540 

3306 

ECT 11-2 X + || 

2512 

3270 

CEM-LMTO 

1820 

2470 

CEM-EMP 

1660 

2310 

MD/MC-CEM-LMTO 

2170 

3010 

MD/MC-CEM-EMP 

1960 

2810 

EAM 

1544 

2082 


as a measure of the smoothness of the surface, 
then the unrelaxed (210) surface is smoother than 
one would assume from its roughness value (B r ). 
Table 8 shows the borocity as a function of planes 
included for some low-index fee faces. 

The lower borocity of the (210) face accounts 
for the fact that, within the hard spheres scheme, 
the (110) face shows a certain degree of trans- 
parency as opposed to the complete coverage in 
(210) faces after a few planes are included. The 
borocity values for the rigid and relaxed Al(210) 
face show some interesting trends. Table 9 dis- 
plays these results. While all ‘solutions’ predict a 
lower borocity than the one corresponding to the 
unrelaxed case, it is rather surprising to see that 
there is little change in borocity between the 
ECT and the experimental values. Although the 
relaxed distribution in each case is different, the 
overall surface effect is quite similar in both 
cases. The comparison between the ECT results 
with perpendicular relaxation only and the fully 
relaxed ones is consistent with the magnitudes of 
the parallel relaxations listed in Table 7. 


Table 12 


Percentage change in interlayer spacing (Sd tJ ) and registry 
shift (A a u ) of Cu(210) and Ni(210) 


±11 

Cu(210) 
ECT II-4 

Ni(210) 
ECT 11-4 

Ad t 2 

- 4.5 ±5.1 

-4.5 ±5.0 

a<* 2>3 

-4.9 ±4.2 

— 5.0 ± 4.1 

Arf,. 4 

±1.0±4.4 

+ 1.1 ±4.4 

Ad 4 . 5 

-2.0 ±5.8 

— 2.0±5.8 

Arf 5 .6 

+ 3.2 ±6.8 

+ 3.2 ±6.7 

Afll.2 

+ 0.0 ±2.4 

+ 0.0 ±2.4 

A a 2,3 

+ 0.0 ±2.5 

+ 0.0 ±2.6 

A a 3,4 

+ 0.5 ±2.8 

+ 0.5 ±2.8 

a<j 4 ,j 

+ 0.1 ±3.5 

-0.1 ±3.4 

A«5,6 

-0.5 ±4.3 

— 0.4 ± 4.1 


Since there is considerable experimental inter- 
est in stepped and kinked surfaces, we include 
predictions of multilayer relaxations for other fee 
metals, Cu and Ni, in order to provide theoretical 
results for future comparison. Table 10 indicates 
the surface energy of Cu(210) and Ni(210) as 
obtained with different approaches, using the 
same notation as in Table 2. 

As with Al(210), we single out the different 
contributions to the surface energy, as computed 
with ECT, in Table 11. Finally, Table 12 displays 
our predictions for the perpendicular and parallel 
relaxations for Cu(210) and Ni(210) usintg ECT 
II-4. 


5. Conclusions 

Based on the analysis of the experimental and 
theoretical results for multilayer relaxation of the 


Table 11 


Surface energy of Cu(210) and Ni(210) (a) (in erg/cm 2 ) and the 
formulations of the bond-compression effect (see text) 

different ECT 

contributions as 

obtained with different 

Surface 

ECT 

(T 

<j | 


0*3 

04 

Cu(210) 

ECT II-4 

2540.16 

2508.93 

29.30 

0.96 

0.96 


ECT 11-2 

2512.25 

2465.42 

43.09 

1.76 

1.97 


ECT I X +|| 

2559.19 

2510.25 

46.98 

1.00 

0.95 

Ni(210) 

ECT 11-4 

3305.79 

3264.75 

37.78 

1.82 

1.44 


ECTII-2 

3269.73 

3208.40 

55.11 

3.29 

2.93 


ECT I ± +|| 

3329.90 

3265.54 

60.91 

1.97 

1.48 
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Al(210) surface, we addressed several issues re- 
garding the implementation of semiempirical ap- 
proaches to the study of such phenomena. We 
noted that surface relaxations other than recon- 
struction involve small energy changes which may 
be difficult to determine accurately considering 
approximations used and the precision of the 
input parameters. In examining these issues for 
ECT applied to the Al(210) surface, we find that 
we obtain different results for experiment and 
considerable uncertainty in the theoretical pre- 
dictions. The quality of the equivalent crystal 
theory results facilitates the discussion on the 
influence of several factors, both internal and 
external, on the ensuing results: the quality of the 
experimental input used, the mechanisms present 
in the algorithm for describing the behavior of 
the system and the analysis of the results in terms 
of relevant properties associated with the system 
under study. Therefore, at present we feel that 
conclusions based on theoretical values of surface 
relaxations are, at best, only meaningful to the 
extent that they refer to relaxation patterns and 
relative magnitudes. 
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